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Abstract. The gravitational instability of expanding shells is discussed. 
Linear and nonlinear terms are included in an analytical solution in the 
static and homogeneous medium. We discuss the interaction of modes and 
give the time needed for fragmentation. Masses of individual fragments 
are also estimated and their formation rates and the initial mass function 
are derived. Results of simulations are compared to observation. 
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Introduction 



' The time evolution of the gravitational collapse of perturbations in decelerat- 

ing, isothermal shocked layers has been examined numerically and analytically 
using linearized equations by Elmegreen (1994). In this paper, we continue with 
nonlinear analysis of the gravitational instability of spherically symmetric shells 
Q-f expanding into the stationary homogeneous medium. We modify the approach 

adopted by Fuchs (1996) who described the fragmentation of uniformly rotating 
self-gravitating disks. The inclusion of higher order terms helps to determine 
with better accuracy when, where and how quickly the fragmentation happens. 
Formation rates and the mass distribution function of fragments are derived 
from simulations and compared to the observed mass spectrum of molecular 
clouds. 



2. Hydrodynamical and Poisson equations on the surface of a thin 
shell 

Equations as derived by Wiinsch & Palous (2001) are 

dm 

— + (X7,mv) = AV Po , (1) 



1 d(mvj 2 



A dt 



c 2 VS - £V$, (2) 



A$ = 4vrGS5(z) , (3) 

where V is the expansion velocity of the shell, m is mass in the area A on the 
surface of the shell, and v denotes a two dimensional velocity of surface flows. 
S is the surface density, c is the constant isothermal sound speed inside the cold 
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Figure 1. Left: the evolution of the maximum surface density per- 
turbation in the linear and nonlinear case. Right: the mass spectrum 
of fragments. The straight line is the power law fit of the decreasing 
part of the spectrum mj^. 



shell, $ is the gravitational potential generated by the mass distribution in the 
shell, G is the gravitational constant and S(z) is a delta function of the space 
coordinate z perpendicular to the surface of the shell. 

2.1. The linear solution 

The solution of linearized Eq. (1-3) may be written in the form of exponential 
functions for the perturbations: Si_£ (perturbed surface density) and v ~ e lujt . 
The shell starts to be gravitationally unstable when 




is purely imaginary and negative. R denotes the radius and So the unperturbed 
surface density of the shell, to is the angular frequency, in the case of instability 
iuj gives the perturbation growth rate, and n the wave number of perturbations. 
The time evolution of So and of Si_l, is shown in Fig. 1 for shell expanding in 
the homogeneous medium of the constant volume density 1 cm~ 3 and average 
atomic weight 1.3. The total energy released is 10 53 erg and c = 1 km s _1 . 

2.2. The nonlinear analysis 

Wiinsch & Palous (2001) include quadratic terms into equations for the evolution 
of perturbations. The time evolution of the perturbed maximum surface density 
including nonlinear terms, Si_jv£, is given in Fig. 1. As another result of the 
nonlinear analysis Fig. 2 shows the spatial distribution of the surface density 
and velocity vectors of the surface flows resulting from mode interactions. 



3. Time evolution of fragments and their mass spectrum 

After the time when the shell becomes unstable, the mass of the fragment de- 
creases as its size decreases. Later, the inflow of mass begins to dominate (surface 
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Figure 2. Left: the spatial distribution of the surface density at t = 
50 Myr. Right: velocity vectors of the surface flows at the same time. 

flows are shown in Fig. 2) and the fragment mass grows. Masses of individual 
fragments are between a few times 10 3 M Q and 1.5 x 1O 6 M0 with the highest fre- 
quency around mf rag ~ 10 4 M Q . Another consequence of the dispersion relation 
(4) is that the most massive fragments, corresponding to the smallest unstable 
wave number rj, form very slowly. Therefore, the most massive clouds are cre- 
ated at the latest stages of the shell evolution. The mass spectrum of fragments 
is shown in Fig. 1. Its decreasing part can be approximated with a power law 
dN/ dm f rag ~ m f*rag w here a = —1.4. This is quite close to the observed mass 
spectrum of GMC in the Milky Way. Combes (1991) gives a = -1.5. NAN- 
TEN survey of the CO emission of the LMC (Fukui, 2001) gives steeper slope: 
a = —1.9, which may be connected to higher level of random motions in the 
LMC compared to the Milky Way, which restricts the formation of late time 
massive fragments in the LMC and steepens their mass spectrum. 
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